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The possibility to induce predetermined coherent quantum dynamics by controlling only the dissi¬ 
pative environmental parameters (such as temperature and pressure) is studied using the combined 
optimal control and environment engineering frameworks. As an example, we consider the problem 
of transforming an optically excited donor state into free charge carriers via intermediate higher- 
lying bridge state(s), with a view to solar energy conversion. In this context, vibrational bath 
engineering allows to promote fast, directional charge transfer and to suppress recombinative losses. 
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I. INTRODUCTION 

Understanding the charge and energy transfer in bi¬ 
ological and synthetic molecular systems is essential for 
the development of efficient ecological energy sources and 
molecular electronics nano-devices Hi- The dynamics 
of these systems belongs to a gray area between quantum 
and classical physics, where the efficacy and sustainabil¬ 
ity of the quantum coherent mechanisms is challenged by 
strong couplings with complex environments. Determin¬ 
ing the optimal composition of coherent and incoherent 
couplings is challenging theoretically and experimentally. 
In particular, the role of the quantum coherences in fast 
and lossless charge transfer (CT) in biological photosyn¬ 
thetic complexes Q remains the subject of debate 0- 

Tailoring the charge transfer system and its environ¬ 
ment to be the most efficient in converting the inci¬ 
dent incoherent solar radiation into free charge carriers 
is of great fundamental and practical interest. However, 
this problem lies outside the scope of the conventional 
quantum optimal control (OC) theory, which is primar¬ 
ily focused on guiding the dynamics of predefined sys¬ 
tems by sha ping the external electromagnetic impacts 
(see e.g. [H, Il4j and references therein). At the same 
time, new approaches to engineer the environment in dis¬ 
sipative dynamical systems for quantum state prepara¬ 
tion and dynamical control were recently suggested [H- 
and experimentally demonstrated [ 23 . In particular, 
they can be efficiently used to prepare highly entangled 
quantum states with high accuracy [ITI - I^ as well as to 
rapidly switch on and off multiple decay channels with 
ultrafast precision [ 2 ^. These studies together manifest 
that the external synergetic synthesis of coherent con¬ 
trols out of disorganized external energy inflow (as repre¬ 
sented by laser generation) can be effectively substituted 
with the local self-organization process in the immedi¬ 
ate quantum system’s environment. Moreover, the envi¬ 
ronmental controls potentially allow to steer the system 
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into regions of the Hilbert space which are out of reach 
for optical coherent control . These new achieve¬ 

ments suggest considering the problem of effective charge 
transfer system design as a problem of integral quantum 
optimization of both junction properties (“Hamiltonian 
structure controls” ISOj ) and reservoir parameters (“in¬ 
coherent controls” [ 23 ). 

In this paper we show that even relatively simple en¬ 
vironment engineering allows to get remarkable bene¬ 
fit from constructive cooperation of resonant transitions 
and decay processes. Specifically, considering the case of 
charge transfer events, we propose a robust and flexible 
mechanism to compensate the unfavorable intrasystem 
electronic arrangement when the charge transfer is orig¬ 
inally strongly suppressed due to unfavorable energetics. 
Our results are based on a simple but very general 3- 
electronic-state model of a junction in which the initial 
excited state of the donor decays into free carriers via an 
intermediate charge-transfer state. Despite its simplicity, 
the 3-state model is useful for describing single-molecule 
junctions molecular bridges between quantum dots or 
molecules and semiconductors singlet-triplet exciton 
fission [3, II] and conversion Q , as well as for semiquali- 
tative analysis of the photosynthesis M- Its extensions, 
including multiple intermediate bridge states, also can 
account for the effect of quantum electronic interferences 

013. 

The paper is organized as follows. First, we present the 
mathematical details of our model and the optimization 
procedure. Then, we qualitatively and quantitatively an¬ 
alyze the outcome of numerical simulations. The con¬ 
cluding part contains the brief summary of the key re¬ 
sults. 


II. THE MODEL DESCRIPTION 

In our model the donor excitation g—^-e by a random 
daylight photon evolves with probabilities Pi and Pu to¬ 
wards two channels (Fig. [T]): (I) electron-hole recombi¬ 
nation back to the ground state g of donor, or (H) the 
two-step charge separation e—Here, oxidation of 
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the donor is followed by slow incoherent decay of a charge 
transfer state ct into a reduced state a of acceptor. Our 
goal is to maximize Pu while suppressing the recombina¬ 
tion. 

All four involved electronic states are assumed to be 
weakly coupled with the modes of a cold Markovian 
phonon bath except for two modes xi and X 2 with fre¬ 
quencies and V 2 - These two modes together with the 
mentioned electronic states form the basis |z, , Ux^) of 

our reduced Hilbert space ^ where z=g, e, ct, a and Ux^ 
are the associated vibration quantum numbers. 

In sequel we consider for simplicity the case of natural 
broadband illumination with constant spectral intensity 
in the relevant frequency range and neglect any sequential 
photon absorption effects which are unlikely on the real¬ 
istic decay timescales in channels (I) and (II). At these 
conditions, the photon absorption at t =0 will lead to ini¬ 
tial population of excited levels |e, Ux ^, proportional 
to square of their overlap integrals with the equilibrium 
state ptd- The resulting excited state 

p|t=ooc ^ |e, n'x^ , (e, rix^ , | x 

J.1 5 ^2 ’ ’ a:2 

Tr [TT[\e,nx^,nx^) {e,n'x^,n' WTrlptd]] (1) 
will further evolve according to the Liouville equation: 

|^=-^[i?,p]+ 2 diss[p]. ( 2 ) 

The Hamilton term of eq. m accounts for the charge 
transfer e-o-ct with coupling constant g and electron- 
phonon interaction with modes Xi and X 2 '- 

H=g|e) (ct|-t- ^ |z) (z| ^ {(dl+afe,z)x 

z=g,e,ct ^ k—1 

(ak+ak,z) + ^} +/3z I +h.c. 

(3) 

Here Bz are electronic states energies, is the phonon 
annihilation operator for the fc-th mode in electronic 
state e, and the constants ak,z=\/ ”^ 2 h^ {xk,z—Xk,e) are 


the scaled dimensionless shifts {xk,z—Xk.e) of the vibra¬ 
tion potential energy minima Xi^z relative to ones in the 
state e. The last term on the right-hand side of eq. © 
describes vibrational intermode coupling with coupling 
constant /3z- 

The dissipation part, Sfdissip], in © has form: 

Sfdiss =rsepSS (Ca, ct) [p] TbrecSS (Csje ) [p] + 

2 

^7fcSe(Cvib.fc)[p]+ ^ ra2(C'z.z)[p], (4) 

k—1 z—g,e,ct 

where S£(C') are conventional Lindblad terms: 

SE{C)[p]=CpC^-^{C^Cp+pC^Cl (5) 

and 

Cz2,zi— 1^2) (^l|i C^vib,fc— ^ ^ |z) (z|(n/c-l-Offc^z)- (6) 

z=g,e,ct 


The first two terms in © represent the charge recombi¬ 
nation e~~>g and charge separation ct~~>a decay channels 
with the corresponding decay rates Tree and Tgep. The 
remaining terms describe the thermal equilibration of the 
vibrational modes caused by weak Markovian collisional 
interaction with the cold bath {kT<^hvk)- They include 
both inelastic contributions with the rates 7 ^ and the 
elastic state-specific exponential transversal decays of the 
vibronic coherences with the rates Tz. 

The dissipation dynamics in © are described in terms 
of zero-order states |z, and do not account for 

the effects of intermode and electronic couplings. The 
reliability of this approximation will be supported below 
by the numerical results. 

We will consider the “worst-case” situation shown in 
Fig. [1] when 1) the charge transfer state ct is above the 
state e: A=Ect—Eg>0; 2 ) electronic coupling is small: 
g<^A; 3) the electron decay rates are equal and slow: 
rsep=ri.ec‘Cp. At these conditions and in absence of vi¬ 
bration modes xi and X 2 the value of Pn is small (at least, 
is evidently less than 50%). Our goal is to increase the 
outcome of channel (H) by adjusting the phonon “dress¬ 
ing” which can be posed as an OC problem: 


Pii 

- > max 

Pi g,VkP^,ak,^P^ 


(7) 


where the parameters g, i^fc, 7fc, Fz, a^^z and /3z 
(z=g, e, ct; A:=l,2) of vibrational environment are 
treated as “controls”. The problem 0 was solved by 
numerical propagation of eq. © in conjunction with 
the genetic algorithm. The details are presented in 
Appendix. 


III. RESULTS 

A. Qualitative description 

Let us first outline the rationale behind the optimal 
“control” parameters found in the simulations. The char- 
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acteristic feature of all converged results (up to exchange 
of the vibration mode indices) is the following relations 

(Fig.EK): 


\/:^\=h{KV 2 -vi) (k = 1,2,3....); (8a) 

7i>0; 72=0 (8b) 

The origin of relations ([8]) for the case of k= 1 is clari¬ 
fied in Fig. HKb). One can see that the quantum of en¬ 
ergy from the mode X 2 (xi) is required to support the 
forward (backward) resonant transition e—>-ct (ct—>e). 
The relations (I8bl) imply that the vibrational energy re¬ 
source for forward transitions is conserved whereas the 
excitations in xi-mode “leak away” into the bath modes. 
Thus, if the geometry of the vibrational potential energy 
surface is such that the initial optical excitation is accom¬ 
panied by strong vibrational excitation of the X 2 -mode, 
then the forward transitions are resonant and fast while 
the backward transfer is off-resonant and slow. Thus, the 
population will be effectively locked in the ct state until 
all the optically induced energy supplies in X 2 -mode will 
be consumed. 


B. Numerical illustration 


As an example, we present the numerical results for 
the case of ^=^0 and rsep=rrec=0.01. In order to get 
a clue into the excited state dynamics it is instructive to 
consider the model case in which the incident light uni¬ 
formly populates the cluster of N 1 XN 2 vibronic states: 


p\t=o— 


1 

N 1 N 2 


Ni SN2+N2 

E E I®; ’^Xi : ^.3,2) (s, Tlxi 1 ^X2 I ■ 

raxi=0 nx2='57V2 


(9) 

All the identified solutions of the OC problem (O with 
initial excited state (O satisfy eqs. dS]), where the value 
of K depends on the choice of summation limits in Q. 
Specifically, the optimal parameters for choice (5A^2=3, 
Afi=ll, and A^2=13 correspond to fc=3: ai^ct=0.5, 



Figure 2. (a) Scheme of optimal conditions for effective uni¬ 
directional electron transfer from the e to ct state (for k= 1). 
(b) Mechanism of population locking in ct state: forward 
e—>^ct resonant transfer is possible whenever there exist vibra¬ 
tional quanta in X 2 mode; backward ct—>-6 resonant transfer 
is locked just after relaxing the x\ mode. 


a2g=-2.69, a2ct=-3.02, j^i=9.28, j^ 2=13.09, f=2.35, 
y,=0.26, re=ret =72 = 0, /lg=0, /3e = -0.11, /Ict=0.17. 

The contributions of the individual levels |e, , na; 2 ) 

to the overall performance Pn ~ 0.85 of this solution 
are detailed in Fig. [3^. One can see that the probabil¬ 
ity of reaching the acceptor state increases abruptly once 
71x2 (he. once the resonant e—>ct transition becomes 
allowed) and slightly decreases with increasing n^i (re¬ 
call that Xi-quanta “promote” the backward transitions 
ct-^>e). We have checked that the result is virtually in¬ 
sensitive to the presence of initial coherences between 
states |e, , 712 , 2 ). the same time, the performance 

degrades substantially in the absence of the intermode 
couplings /3e and /3ct (Fig. [51 d). Comparison with Fig. [3^1 
shows that these couplings enhance the resonant transfer 
rate in the cases when the transitions are suppressed by 
small values of Franck-Condon overlap integrals. 

The efficacy of charge separation is also very high in 
the case of broadband incoherent solar energy absorp¬ 
tion (see eq. ©)• In this case the maximal performance 
Pn — 0.85 can be achieved with ai.g=0 at low tempera¬ 
tures, when the equilibrium vibronic state ptd coincides 
with the ground state |g,0,0). However, the results of 
calculations summarized in Fig. 2] allow to predict sub¬ 
stantial efficacy losses at higher temperatures when the 
thermal excitation of mode X 2 is non-negligible. 


Fig. [ 5 ] provides a complimentary way to elucidate the 
charge separation mechanism in terms of the evolution 
of the individual eigenstates of the Hamiltonian ©. 
One can see that the major “feeders” of channel (I) are 
strongly localized on the donor. More detailed analysis 
shows that all of them are composed of such zero-order 
states Ie, 712 : 1 , 713 , 2 ) fFat nx^<K,, so that the forward trans¬ 
fer e—>-ct is energetically forbidden. In contrast, the vast 
majority of the remaining excited states predominantly 
evolve towards the acceptor side regardless of their struc¬ 
ture and localization. Due to thermal relaxation, the 
populations of these states at times — — will 

”^2 i sep,rel 

be accumulated in the lowest energy “drain” state of the 
upper group (indicated by an arrow in Fig. [5]). This 
state has pronounced charge transfer character and is es¬ 
sentially decoupled from the donor side. For this reason, 
its decay dynamics can be accurately represented within 
the zero-order basis approximation. These observations 
confirm the applicability of the model assumptions made 
above. 

We remark, interestingly, that the “drain”-type de¬ 
coupled eigenstates that sustain the radiationless decay 
within the strongly mixed vibronic systems, have been 
discovered in realistic molecular models and studied in 
context of optimal control [^ . Their remarkable prop¬ 
erties have been related to the phenomenon of scars of 
periodic orbits [3^ . 
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Figure 3. (a) The probability of various excited states |e, to evolve towards the channel (II) for optimized system 

parameters (see the text), (b) The same dependence with disabled vibrational couplings: /3e=/3ct=0. The nodes of red line 
show the relative absolute values of the overlap Franck-Condon integrals for |e, , n^j) —>■ |ct, —k) transition. 


Pu 



Figure 4. Overall efficacy of the light-induced charge separa¬ 
tion as a function of the vibrational state \?,,nxx,nx 2 ) of the 
parent state. The excitation is considered within the broad 
spectrum approximation. Only the dependence on the 
vibrational quantum number is shown for fixed na;^=0; the 
dependence on is negligible. The “controls” were opti¬ 
mized for launching from the ground vibrational state of g 
with na;j=na; 2 = 0 . 



Eigenstate energy 


Figure 5. Probability of reaching the outcome (II) starting 
from the different eigenstates of the Hamiltonian (|3]). The 
localization of each eigenstate is indicated by color codes as 
clarified in the inset. 


IV. SUMMARY AND CONCLUSION 

In this work we proposed an approach to (incoherently) 
controlling daylight-induced charge separation, hence a 
potential route to solar energy conversion. The control 
mechanism is 

• universal. It can be applied without any substan¬ 
tial change to stimulate a variety of other types 
of unidirectional quantum dynamics (e.g. energy 
transfer, proton transfer etc.) in different donor- 
bridge-acceptor systems; 

• effective, irreversible and fast. In fact, it is capable 
to enforce even counterintuitive sustainable popu¬ 
lation flow from the lower-energy to higher-energy 
electronic state; 

• essentially insensitive to the spectral characteris¬ 


tics of the incident light and hence potentially well- 
suited to harvesting broadband solar energy. 

The substitution of lasers with readily available cheap 
energy sources is a fascinating challenge with obvious 
ramifications. Based on the results of this work we pro¬ 
pose the following concepts for achieving this goal: 

• The concept of a vibrational reservoir as a dual-use 
resource to achieve unidirectional quantum dynam¬ 
ics: The associated coherent mechanisms promote 
the resonant forward e—>ct transfer while the in¬ 
coherent decay channels are responsible for pop¬ 
ulation locking in the target state ct. From the 
energy standpoint, this concept exploits the fact 
that optical excitation, in general, brings both the 
electronic and the vibrational subsystems into non¬ 
equilibrium, and subsequent evolution of the latter 
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imposes coherent feedback on the former. As a re¬ 
sult, all the essential quantum coherences are cre¬ 
ated locally without use of external coherent energy 
resource. 

• The concept of marrying quantum optimal con¬ 
trol methods with environment engineering to de¬ 
termine the characteristics of the reservoir which 
will lead to a desired outcome. In other words, the 
quantum optimal control is put in the framework of 
the traditional chemistry methods, so that not the 
laser fields but the temperature, pressure, parame¬ 
ters of solvent, chemical composition etc. become 
the key quantum control parameters. 

One of the interesting questions that this study opens 
is whether the proposed mechanism is also exploited by 
natural light-harvesting complexes. Specifically, it pro¬ 
vides the natural framework for converting any initial 
coherences which can be encoded in the certain vibra¬ 
tional modes by laser excitation into electronic coher¬ 
ences similar to observed in the recent 2D-spectroscopic 
experiments ii- The analysis and verification of this 
possibility would carry us beyond the scope of the present 
paper. 
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APPENDIX: THE OPTIMAL CONTROL OF 
CHARGE TRANSFER 


1. Numerical estimation of the probabilities Pi and 

PlI 

The definition of Pi and Pu implies that these quanti¬ 
ties are feasible physical parameters only until the elec¬ 
tron dynamics remains confined in the Hilbert space 
"X. In practice, this corresponds to timescales T of or¬ 
der T'^min (Tgep, Prec)”^- The key simplifications (one- 
electron picture, absence of sequential absorption and re¬ 
combination a—s-g, Markovian phonon bath approxima¬ 
tion) embedded into the Liouville equation ([2]) are also 
justified for this time window. However, the character 
of these simplifications ensures that the populations of 
states g and a formed on the timescales of T will re¬ 


main “frozen” during further “unphysical” propagation 
of eq. (HI). In particular, this means that: 

p|t=oo=Ti6'i+Tii£'ii, (10) 

where gi= |g, 0,0) (g, 0,0| and Qii= |a, 0, 0) (a, 0,0| are the 
linearly independent steady eigenstates of the total Li- 
ouvillian S£[©]=^[H, 0]-|-SCdiss©: SS[gpii]=0. Denoteas 
'I'l and the associated normalized right eigenvectors: 
vl/ise[o]=0; Tr[p^vl- 5 ]= 1 ; Tr[p^^4VI-5]=0 (^,x=I,H). The 
last three relations lead to the following identity: 

Tr[^ii (Se 0 +PII Tr[6©] + qi Tr[©])]=Tr[0©], (11) 

where O is an arbitrary Hermitian operator such that 
the superoperator in the round brackets is full-rank (e.g. 
0=Qii). Multiplying (ITOl) by ^'n and taking trace gives: 

Pii=Tr[^'iip| ]=Tr[^'iip|to]; Pi=l-Pii, (12) 

where the fact of time-independence of is used. 

Eqs. (fTTl) and (TT^ constitute the complete system for 
determining Pi and Pu. In simulations we represented 
eq. dTT]) in appropriate operator basis as a system of lin¬ 
ear equations Lij'$j=Oi relative to and then solved 
it using the conjugate gradient method with a diagonal 
Jacobi preconditioner M=diag{Li,i, ^ 2 , 2 , ■•■} 

2. Solution of the extremal problem Pn/Pi—max 

The approximate extrema of the objective functional 
^max were searched using the genetic algorithm. In 
our calculations we used the genotype which includes 11 
parameters: 5 , 1 / 1 , 1 / 2 , 71 , 72 , Te, Tct, afc,e, afc,ct, /3e and 
Pet as well as its extended 13-parameter variant where 
the vibrational frequencies vi and V 2 are allowed to be 
the functions of electronic state z. We found, however, 
that the maximal efficacy is achieved in absence of de¬ 
pendences j/i_ 2 (z). The optimization was performed on 
the population of 25 individuals. Each subsequent gen¬ 
eration includes the best 34% of parents, another 33% 
are obtained via random mutation of the 20% of genes in 
seeding individuals and the rest result from crossover of 
20% of genes between two seeding individuals. The seed¬ 
ing individuals were randomly selected among the best 
50% of parents merged with extra 3-4 randomly gener¬ 
ated candidate solutions. 

After reaching substantial convergence the best solu¬ 
tion was further refined using the standard gradient de¬ 
scent method. 
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